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Abstract 

Pulsars are rapidly rotating neutron stars and are the outcome of the 
collapse of the core of a massive star with a mass of the order of or larger 
than eight solar masses. This process releases a huge gravitational energy 
of about 10 53 erg, mainly in the form of neutrinos. During the collapse 
the density increases, and so does the magnetic field due to the trapping 
of the flux lines of the progenitor star by the high conductivity plasma. 
When the density reaches a value of around 10 12 gcm~ 3 neutrinos become 
trapped within the protoneutron star and a neutrinosphere, characterized 
inside by a diffusive transport of neutrinos and outside by a free stream¬ 
ing of neutrinos, is formed and lasts for a few seconds. Here we focus on 
the structure of the neutrinosphere, the resonant flavor conversion that 
can happen in its interior, and the neutrino flux anisotropies induced by 
this phenomena in the presence of a strong magnetic field. We present 
a detailed discussion in the context of the spherical Eddington model, 
which provides a simple but reasonable description of a static neutrino at¬ 
mosphere, locally homogenous and isotropic. Energy and momentum are 
transported by neutrinos and antineutrinos flowing through an ideal gas of 
nonrelativistic, nondegenerate nucleons and relativistic, degenerate elec¬ 
trons and positrons. We examine the details of the asymmetric neutrino 
emission driven by active-sterile neutrino oscillations in the magnetized 
protoneutron star, and the possibility for this mechanism to explain the 
intrinsic large velocities of pulsars respect to nearby stars and associated 
supernova remnants. 


1 


1 Introduction 


The finding of pulsars in the sixties was one of the most astonishing discoveries 
in astrophysics]]]. Soon after their observation they were identified with neu¬ 
tron stars | 21 , whose existence had been theoretically predicted decades before. 
Neutron stars are the densest, most rapidly rotating, and most strongly mag¬ 
netized objects in the Galaxy, features well understood by the standard physics 
of gravitational and elementary particle physics [ 3 ] • They are unique laborato¬ 
ries, which give us access to extreme conditions that are virtually impossible to 
obtain on Earth. Neutron stars are of particular interest for neutrino physics 
because of the key role that these particles play in the very first stages of their 
formationj?]. 

There is a remarkable characteristic of pulsars not satisfactorily understood: 
their large drift velocities with respect to nearby stars. The average value of 
pulsar velocities is in the range of 200 — 500kins -1 , an order of magnitude 
larger than the mean velocity of ordinary stars in our Galaxy (~ 30 km s -1 ) {5]. 
In addition, the distribution of these velocities is not Gaussian. It seems to 
have a bimodal structurejS], with a significant fraction (~ 15%) of the pulsar 
population having velocities higher than lOOOkms -1 and up to a maximum of 
1600 kms -1 . 

Neutron stars are formed in type-II supernova explosions and, given the 
enormous energy liberated during these processes, it is quite natural to look for 
an explanation in terms of an impulse or kick received during their birth. This 
hypothesis is supported by recent analysis of individual pulsar motions and of 
associations between supernovae remnants and pulsars. It is also consistent with 
observed characteristics of binary systems with one or both constituents being 
a neutron star]?]. Despite this evidence, the physical origin of the kick is one of 
the most important unsolved puzzles in supernovae research. If the distribution 
is bimodal, then more than one mechanism could be responsible for the natal 
kicks. Two classes of mechanism have been considered|S]. One of them invokes 
an asymmetric mass ejection during the supernova explosion as the result of 
hydrodynamic perturbations caused by global density inhomogeneities prior to 
the core collapse]^] or convective instabilities in the neutrino-heated layer behind 
the shock [Til. Recoil speeds as high as ~ 500kms -1 have been obtained with 
two-dimensional simulations. However, a recent three-dimensional calculation 
shows that even the most extreme asymmetric collapse does not produce final 
neutron star velocities above 200kms -1 in, rendering the situation unclear. 

Only a small fraction of the gravitational binding energy (E ~ 10 53 erg) 
liberated during a supernova explosion is required to account for the huge elec¬ 
tromagnetic luminosity observed and the ejection of the envelope. More than 
99% of the energy is released in the form of neutrinos, which are copiously pro¬ 
duced when the iron core turns into a neutron core. A small asymmetry (~ 1%) 
in the momentum taken away by neutrinos should be enough to give the nascent 
pulsar an acceleration consistent with the measured velocities. The asymmetric 
neutrino emission induced by strong magnetic fields is the basis of the second 
class of kick mechanisms, and effects such as parity violation 1121 . asymmetric 
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field distributionsfI3|. and dark spots[2I have been considered in this context. 

The density in the interior of a protoneutron star is so high, of the or¬ 
der of the nuclear density, that even neutrinos remain trapped, being emitted 
from an approximately well defined surface, the neutrinosphere. In a magne¬ 
tized stellar plasma there are anisotropic contributions to the neutrino refrac¬ 
tion index |T5I ITC] and, as a consequence, the resonant flavor transformations 
of neutrinos diffusing in different directions respect to the magnetic field occur 
at different depths within the protostar. This phenomenon can induce a mo¬ 
mentum flux asymmetry and is the basis of the neutrino-driven kick mechanism 
proposed by Kusenko and Segre (KS)J7„ [HI El- With standard (active) neu¬ 
trinos the mechanism works when the resonance region lies between the v e and 
i/p (or i/ T ) neutrinospheres, but requires an exceedingly high square mass differ¬ 
ence, in conflict with the existing limits. A variant avoiding such a limitation 
can be implemented in terms of active-sterile neutrino oscillation[20j EH] !221 ■ 

The measurement of the Z-boson width at LEP, limiting to three the num¬ 
ber of light neutrinos having weak interactions, does not preclude the existence 
of additional neutral fermions that do not fill the standard model gauge in¬ 
teractions. They have been invoked to explain several questions in astronomy 
and cosmology j23( and may becomes a necessity if the claimed LSND evidence 
for oscillations were confirmed. The range of the oscillation parameters (mass 
and mixing angle) required by the pulsar kick mechanism overlaps with the 
allowed region for sterile neutrinos being candidates to cosmological dark mat¬ 
ter. Alternative kick mechanisms driven by neutrino oscillations have been 
proposed invoking the possible existence of transition magnetic moments j24). 
non orthonormality of the flavor neutrinos|25j. violations of the equivalence 
principle (2Ej, or off-resonant emission of sterile neutrinos ET.. 

In this article we review the explanation of the pulsar kick in terms of the 
resonant neutrino oscillations. In Section II, we outline the creation of a pro¬ 
toneutron star as a result of a supernova event, and the formation in its interior 
of a gas of trapped neutrinos. The momentum flux carried by neutrinos is an¬ 
alyzed in Section III. Section IV is devoted to the discussion of the neutrino 
oscillations in matter incorporating the effect of strong magnetic fields. The 
implication of this phenomenon as source of the asymmetry in the momentum 
flux transported by neutrinos emitted during the protoneutron star cooling is 
examined in Section V. Numerical results are obtained in Section VI within the 
context of the Eddington model, where the kick effect is explicitly evaluated. In 
the final section we summarize the outstanding features of the mechanism. 

2 Stellar collapse, protoneutron stars, and neu¬ 
trinospheres 

Neutron stars are extremely compact objects born in type II supernova explosions |27). 
as the aftermath of the gravitational collapse of a massive star (M > 8M@) when 
it exhausts the sources of thermonuclear energy and its iron core reaches the 
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Chandrasekhar limit, Mch = 1.44 M 0 . The process is triggered by the pho¬ 
todissociation of iron nuclei, 7 + 56 Fe —> 13 4 He + 4n — 124,4MeV, which 
consumes energy reducing the pressure of the degenerated electron gas. When 
the collapse proceeds, pressure support decreases because of the electron cap¬ 
ture by nuclei. Electron neutrinos are copiously produced through the reaction 
p + e —> n + u e and initially they freely escape from the core, but as den¬ 
sity increases the medium becomes less transparent for them. The inner core 
(M « 0.5 — 0.8 Mq) collapses in a homologous way at a subsonic velocity pro¬ 
portional to the radius (v/r = 400—700s _1 ), with the density and temperature 
profiles remaining similar and only the scale changing with time !2S| . The ex¬ 
ternal shells collapse at a supersonic velocity of the order of free fall velocity 
(v oc 1 Isjr). 

During the first part of the collapse the main source of opacity for neutrinos is 
the coherent scattering by heavy nuclei, mediated by the neutral current weak 
interaction. The cross section for this reaction is proportional to A 2 , where 
A is the number of nucleons in the nuclei, and is very effective in trapping 
neutrinos (251. If A„ is the mean free path of neutrinos, then the time for them 
to diffuse out of the core, estimated from a random walk through a distance of 
the order of the core radius R, is tdiff ~ 3f? 2 /A„. On the other hand, the core 
dynamical time scale is of the order of the free fall time t c ~ (Gp)" 1 / 2 . At a 
certain density, of about 10 12 g cm 3 for 10 MeV neutrinos, both times become 
comparable and neutrinos are trapped [30j. /3-equilibrium between e and v e is set 
up and electron neutrinos also become degenerated. In this way, the collapsing 
core is the only place, besides the early universe, where neutrinos are in thermal 
equilibrium. During the collapse the lepton fraction Yl = Y e + Y Ve is kept 
almost constant near the value of the electron fraction at the beginning of the 
trapping (Y e ~ 0.35), with only ~ 1/4 corresponding to neutrinos|3|. 

The degenerate gas of trapped neutrinos diffuses towards the exterior, reach¬ 
ing regions of decreasing density where the collisions are less frequent, until 
they finally decouple from the star. Although this process is continuous, it is 
convenient to introduce an effective emission surface called the neutrinosphere, 
characterized as the surface where the optical depth becomes approximately 
one. More accurately, the radius of the neutrinosphere R. Vf is defined by the 
condition 



( 1 ) 


The medium is considered opaque for neutrinos inside the neutrinosphere and 
transparent outside. We can also define the energy sphere, where the energy 
exchanging reactions freeze out while energy conserving collisions may still be 
important. 

The collapse stops suddenly when, after a fraction of seconds, densities of 
around 3 x 10 14 gcm -3 , exceeding nuclear density, are reached pm. Then, a 
shock wave is built up near the interface of the halted inner core and the su¬ 
personically free falling outer part. The shock wave propagates outwards and 
eventually ejects the stellar mantle of the star. At its passage through the neu- 
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trinosphere nuclei are dissociated decreasing the neutrino opacity and causing a 
recession of the neutrinosphere. The v e produced by the reaction p + e —> n + v e 
can escape freely building up a short (< 10 msec) v e burst, called the “neutron- 
ization burst” or “deleptonization burst”. The star remnant gravitationally de¬ 
couples from the expanding ejecta and a protoneutron star has been formed 1321 . 

The neutrinos continue trapped due to the interaction with the residual 
particles of the medium, mainly nucleons, electrons, and neutrinos. For v e 
the dominant reactions are the (3 and /J-inverse, which mantain them in local 
thermodynamical equilibrium j33| . In this case, the distinction between neutri¬ 
nosphere and energy sphere is not crucial. The electron neutrinos also interact 
with nucleons through the neutral current, but the cross sections for these re¬ 
actions are smaller. Although the energy dependence of the process involving 
v e and v e is the same, in a protoneutron star the number of neutrons is larger 
than the number of protons, and therefore the v e suffer more collisions than the 
v e - As a consequence, the neutrinosphere of the v e is a bit further out than the 
one of the i/ e , and the mean energy per emitted particle is smaller. After the 
achievement of (3 equilibrium, the total entropy is conserved and the collapse 
becomes adiabatic. 

In the case of the nonelectron (active) neutrinos and antineutrinos, the ab¬ 
sence of muons and taus does not allow the establishment of chemical equilib¬ 
rium with the stellar plasma. Therefore, they do not transport leptonic number 
outside the star and their chemical potentials vanish. Their interactions with the 
background are the same and no distinction needs to be made between v lt and 
v T . Deeper in the protostar the main interactions for these neutrinos are pair 
processes such as MM' MM'v^^v^^, e~e + etc, where M denotes 

a nucleon, and the scattering reactions with nucleons and electrons. The scat¬ 
tering on electrons are less frequent, but the energy transferred is greater. This 
process dominates the energy exchange up to the energy sphere where it freezes 
out. Beyond this radius, muon and tau neutrinos still scatter off neutrons until 
they reach the neutrinosphere. The cross section for Mv^ T M tv T is smaller 
than the one of the [3 reaction, and hence the radius of the muon (tau) neu¬ 
trinosphere R v T is shorter than the radius of the electron neutrinosphere R Ue . 
The protoneutron star contracts, cools, and deleptonizes due to the emission of 
neutrinos and antineutrinos of all the flavors in a scale of time of the order of 
10 sec. The released energy is equally shared by all neutrino and antineutrino 
species. This phase is known as the Kelving-Helmholtz phase and ends up in 
the formation of a neutron star. During this time the radius of the protostar 
varies from around 30 to 10 km, while the mass remains of the order of 1.4 Mg. 
After approximately 1 minute the neutrino mean free path becomes comparable 
to the stellar radius and neutrino luminosity decreases quickly. 

Another important property predicted for neutron stars is that they should 
have an extremely strong magnetic field. Since the conductivity of the iron core 
is very large, the magnetic flux is essentially trapped and it is conserved as it 
collapses to form a neutron starj^J. Adopting for the initial magnetic field the 
typical values observed in some white dwarfs, from the flux conservation one can 
easily estimate that the surface field of a protoneutron star has to be in the range 
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of 10 12 — 10 13 G. This is supported by evidence inferred from measurements of 
the spin-down rates and observations of the X-ray spectra features pT5|. Magnetic 
fields of the order of 10 15 G or stronger can be generated by dynamo processes, 
and in the central regions intensities up to 10 18 G have been considered[&>}. 


3 Momentum flux 


In a protoneutron star the momentum flux from the core to the external regions 
is dominated by neutrinos, which provide the most efficient mechanism to re¬ 
lease the energy generated during stellar collapse. As explained in the previous 
section, matter densities are so high that neutrinos are trapped forming a de¬ 
generated gas at quasi-equilibrium diffusing radially. The statistical description 
of this gas is made in terms of a distribution function f V( (x, k, t). If the action 
of external fields is negligible, then the neutrino momentum between two con¬ 
secutive collisions remains constant. The change of the distribution function is 
due to neutrino collisions with the particles in the medium and is described by 
the Boltzmann equation 

^ i =C'(/„), (2) 

where C(f Vt ) is the collision integral and 

tt-tM-*/- (3) 

When the regime is close to equilibrium we can approximate the collision 
integral as follows 

f _ feq 

C{f ve ) — - Vt Jvi , (4) 

where t Vi = is the average time between collisions and 

= i + e (fc-^)/T ( 5 ) 

is the Fermi-Dirac distribution function for the gas of relativistic neutrinos 
(E Ue = k = |k|) with chemical potential /i„ f and temperature T. In the static 
situation ( df„ e /dt = 0), Eq. (J2J) with C(f Ve ) given by Eq. (@J) can be solved 
iteratively. Proceeding in this way we obtain 

( 6 ) 


which constitutes the diffusion approximation for the distribution function. The 
second term in this expression gives the correction to equilibrium and allows the 
transport of the energy and momentum towards the exterior of the protostar. 

Once the distribution function is known the neutrino energy-momentum ten¬ 
sor can be computed as 


T » v = [ 

Vi J (2t r) 3 fc° 


(7) 
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From this formula, for the energy density we have 


Uvi = = (2^)3 / fc3rffcsin ^ r-T 4 , (8) 

where 77 ^ = p, Vl /T is the degeneracy parameter and F :i [rj vi ] is a Fermi integral 

pOO 

F n[rjvi]= / dxx n [l + exp(x - r/^)] _1 . (9) 

Jo 

At each point within the protostar we choose a local frame ( x', y\ z') with r = z', 
such that k = cos O' z' + sin O' cos <f)'3c' + sin O' sin <p'y'. 

The stress tensor can be expressed in terms of the energy density as 

_ §ij _t!L _ ( 10 ) 

3 

The equilibrium component of the distribution function is isotropic and gives 
a null contribution to the momentum flux ( F ’* = Tj;®), which is generated 
exclusively by the diffusive term: 

1 r i r°° 

F„, = ^3 J k 2 dk sin O'dO'dtf k f Vl = jf k 3 \„ e \7f%dk . ( 11 ) 

To complete the calculation of F Vl we need the mean free paths of the neutrinos 
or equivalently, the opacities k v , . 

The neutral current reactions equally affect all flavors. Since the energies 
available in the system are much smaller than the mass of the r lepton, the v T 
does not participate in the charge-current interactions. For the same reason, 
and for simplicity, we also ignore the presence of muons. Accordingly, the main 


sources of neutrino opacities are 

vt + n ^ vt + n , (12) 

ve+p^ve+P, (13) 

v e +n^e~ +p. (14) 

In terms of the cross sections for these reactions the A” 1 are given by 

K, = Kt + Merabs, (15) 

K^ t = N n<?n + NpCTp , (16) 

where N ntP are the neutron and proton number densities and 

^ ^ (! + 3^)fc 2 , (17) 

o p = a n 1 - 8 ^ (1-2 sin 2 0 W ) , (18) 

1 + dg A 

<7 abs = 4:<J n . (19) 
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Here, sin 2 Ow — 0.23 and gA — 1-26 is the renormalization for the axial-vector 
current of the nucleons. In Eq. m we neglected contributions of order ( m n — 
m p )/k 1 where m n and m p are the neutron and proton mass, respectively. 

From the previous formulae we immediately see that 

— x u ^k p , ( 20 ) 


where the coefficients x Vf are constants. By making m n = m p and taking 
the typical values Y n « 0.9 and Y p ss 0.1 for the nucleon fractions Y n , p = 
N n , p /(N n + N p ), we get 

x Ue = 3 x lCT 25 MeV“ 5 , (21) 

}T = 0.7 x 10" 25 MeV -5 . (22) 

When Eq. H3) is substituted in the integrand of Eq. o. the momentum flux 
transported by neutrinos diffusing from the core to the exterior of the protostar 
becomes 

P » = lr2 >' <23) 

In what follows we assume that the neutrino chemical potentials are negligible, 
and using Fi[0] = 7r 2 /12 and F 3 [0] = 77r 4 /120, from Eqs. © and J23), we obtain 


U VI 

Y v , 


7tt 2 

240 


rp 4 

1 


72 x vi p 


■ VT 2 


(24) 

(25) 


The last results will be used in Section 6, when we calculate the asymmetric 
neutrino emission in the context of a specific model for the protoneutron star. 


4 Neutrino oscillations in a magnetized medium 

A notable consequence of neutrino mixing is the phenomena of neutrino oscilla¬ 
tions established five decades ago by Pontecorvo pT7113%| . It corresponds to the 
periodic variation in the flavor content of a neutrino beam as it evolves from 
the production point. In the last years compelling evidence of neutrino oscilla¬ 
tions has been obtained in experiments with atmospheric, solar, accelerator, and 
reactor neutrinos pM 14(1) . The most plausible interpretation for these new phe¬ 
nomena is that neutrinos are massive particles and the weak-interaction states 
\vt) (£ = e, p, t) are linear combinations of the mass eigenstates | vj) (j = 1, 2, 3), 
with masses rrij: 

H=U ej \uj), (26) 

where U(j are the elements of a unitary mixing matrix, the analogue of the 
CKM matrix in the quark sector. Consider a beam of v a created at a certain 
point with a 3-momentum k. When the state given in Eq. EU evolves, each 
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of the mass eigenstates acquires a different phase e with Ej = 


+ k 2 


(k = |k|). As a consequence, at a certain distance r = t from the source the 
neutrino beam will not correspond to pure vi but turns partially into other 
flavors. 

The pattern of neutrino oscillations is modified when neutrinos propagate 
through a material medium. The basic reason is that the energy-momentum re¬ 
lation of an active neutrino is affected by its coherent interactions with the par¬ 
ticles that constitute the medium. These modifications are described in terms of 
a potential energy Vi, which is related to the real part of the refraction index ne 
according to Renj = 1—Ve/k and can be calculated from the background contri¬ 
butions to the neutrino self-energy |41|. If, like in normal matter, electrons but 
not heavier leptons are present in the background, then only v e scatter through 
charge current interactions, while all the flavors interact via the neutral cur¬ 
rent. Therefore, the refraction index differs for neutrinos of different flavors and 
they will acquire distinct phases when propagating with the same momentum 
through the medium. This fact can have impressive consequences in the case 
of mixed neutrinos and, under favorable conditions, flavor transformations are 
significantly enhanced in a medium with a varying density even when mixing in 
vacuum is small. This is the essence of the MSW mechanism j32I EH whose large 
mixing realization provides the solution to the solar neutrino problem [¥0122] ■ 

Besides the modifications of the dispersion relation, neutrinos acquire in 
matter an effective coupling to the electromagnetic field by means of their weak 
interaction with the background particles|15|. As a result, in the presence of 
an external magnetic field there are additional contributions to the neutrino 
refraction index and flavor transformations are affected, but in a way that pre¬ 
serves miira,lit,v[44|. contrary to what happens for neutrino oscillations driven by 
transition magnetic moments |4R|. In what follows, we will discuss neutrino os¬ 
cillations in a magnetized medium for the simplest situation of mixing between 
two neutrino species, that is, and another active neutrino Vfj or a hypothetical 
sterile neutrino v s . In any case we can write 



(27) 


where 9 is the vacuum mixing angle, a parameter to be fixed from the experi¬ 
mental results. 

The Dirac wave function for relativistic (mixed) neutrinos with momentum 
k propagating in matter in the presence of a static uniform magnetic field is 
well approximated as follows 25 : 



(28) 


where we use the Weyl representation, is the Pauli spinor of negative helicity 
(a.it (f>- = —</>_) and \(t) is a flavor-space spinor that satisfies 


l d £=M BX . 

at 


(29) 
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The Hamiltonian Eg is a 2 x 2 matrix that, in the basis of the flavor states, 
takes the form 

M 2 

M B =k+ — +Y, (30) 

with 

M = u( m r } ° W (31) 

\ 0 m 2 J 

and ¥ = diag(V), Ve>, a ). 

The effective potential for active neutrinos can be written asjTS:,, H7| 


Ve = be - ce e k • B , 


(32) 


where e > 0 is the proton charge. The coefficients be and ce depend on the 
properties of the thermal background. To order Gf the expression for be is 
independent of B, but this is not true in general for eg. In our analysis we will 
work within the weak-field limit, i.e., B -C /i 2 /2e m, and ignore any possible 
dependence of this coefficient on the magnetic field. The effect of strong mag¬ 
netic fields on the neutrino propagation in a medium have been considered by 
some authors[32]■ For (active) antineutrinos the potentials change in sign, i.e., 
Ve = — Ve , while in the case of sterile neutrinos, since they do not interact with 
the background, we have 

V a =V a =0. (33) 

Eq. 12HI) is the extension of the MSW equation to the situation we are con¬ 
sidering, and was first derived in Ref. [2]. The components Xeif) and Xn,a{t) 
of x(t) give the amplitude to find the neutrino in the corresponding flavor state. 

For a uniform background %(t) = ^ (<r| x(0))fte* a,<t , where w,;(p) are the ener- 

i=l,2 

gies of the neutrino modes in the medium and the vectors ?i i2 are the solutions 
of the eigenvalue condition Egq = cjjC;. For oscillations the relevant quantity 
is 

w 2 — wi = (Am 2 cos 2 9 — 2 fcV ) 2 + (Am 2 sin 20 ) 2 , (34) 

where V(r) = Ve(r) — Vei tS (r) and A to 2 = m\ — m\. When the medium is 
inhomogenous, along the neutrino path (r = t) the elements of V are functions 
of t, and x(t) has to be determined by solving Eq. (T^l with iga time dependent 
matrix. At each point r, can be diagonalized by the unitary transformation 


f COS 9m sin 0 m \ 

V -sin 9m cos 9 m ) ’ 


with the mixing angle in matter 9 m (r ) given by 


tan 29 m 


Am 2 sin 29 
Am 2 cos 29 — 2kV 


(35) 


(36) 


A similar expression is valid in the case of antineutrinos with Am 2 cos 29 + 
2fcV in the denominator. In general V and consequently the mixing angle in 
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matter depends on the magnetic field. As a function of V, sin 2 29 m exhibits the 
characteristic form of a Breit-Wigner resonance. For neutrinos (antineutrinos) 
with momentum k the resonance condition (sin 2 2 9 m = 1) is 


V(r*)=± 


Am 2 

2k 


cos 29, 


(37) 


where r R = r H (k) and the sign + (—) corresponds to neutrinos (antineutrinos). 
Notice that for active-active oscillations the neutral current contributions to 
V = Vi — Vy cancel since they are the same for both flavors. This is not true in 
the case of active-sterile oscillations where V = Ve, and such contributions play 
a relevant role in the resonant transformations for up <—> i/ s (or up <—► i/ s ).The 
efficiency of the flavor transformation depends on the adiabaticity of the process. 
The conversion will be adiabatic when the oscillation length at the resonance is 
smaller than the width of the enhancement region. This imposes the condition 


dV 

dr 


< 



sin 29 


(38) 


The dispersion relation of neutrinos in a magnetized plasma made of elec¬ 
trons, protons, and nucleons have been calculated in Ref. GEj, including the 
contributions due to the anomalous magnetic moment coupling of the nucleons 
to the photon. From the results given there we have 


b e = V2 G f 


Vt = V2 G t 


C e = 2 V 2 Gp 

c ll)T = 2\/2Gf 


N e + - - 2 sin 2 9 W (N p - N e ) - -N, 


+ b 


+ 


- - 2 sin 2 9 W (N p - N e ) - -N n 


g A (l + 2 m p ^j C p - gAZm.n^Cn - C e 
g A (l + 2 m p ^j C p - g A 2m n ^-C n + C e 


(39) 

(40) 

(41) 

(42) 


where n n ^ p are the anomalous part of the nucleon magnetic moments given by 


K n = —1.91- 


2 rrir, 


Kp — 1.79 , 

2m n 


(43) 


and 


bi = V2 Gf 


N Vl — N„ t + 



d 3 p 1 d 
(2tt) 3 2 E dE Jb 


(b = e,n,p). 


(44) 

(45) 


In these equations N e and N VI (Np e ) are the number densities of electrons and 
neutrinos (antineutrinos) of flavor £, respectively and f\, are the distribution 
functions of the electrons, protons, and neutrons. 


11 




















We model the atmosphere of a protoneutron star as a neutrino gas dif¬ 
fusing through a plasma made of relativistic degenerate electrons and classi¬ 
cal nonrelativistic nucleons. In such conditions C e = — /z e / 87 r 2 and C n>p = 

—W r /?>Tm,„ .„ |TC1PTK] . where T is the background temperature and fj, e = (37r 2 fV e ) 
is the electron chemical potential. The neutrino fractions are much smaller than 
the electron fraction Y e = N e /(N n + N p ) ~ 0.1 and, in the context of our anal¬ 
ysis, the contributions to b ue coming from the neutrino-neutrino interactions 
(denoted by b e ) can be neglected. With this approximation in mind, putting 
m p — m n , and taking into account that by electrical neutrality N e — N p , the 
coefficients in the effective potentials become: 


b K ~ 


^/i,T — 


c e ~ - 


Cll.t — 


^ (1 - 3 Y e ) , 

y/2 y ^ m n 

(46) 


(47) 


(48) 


(49) 


It should be noticed that in Ref. m Ce vanishes because all the components of 
the stellar plasma are assumed to be degenerated, and the contributions due to 
the anomalous magnetic moment of the nucleons are not included. 

We now use the effective potentials Vt with be and ce given by Eqs. m-m 
to examine the magnetic field effect on the resonant neutrino oscillations within 
a protoneutron star. In a linear approximation, we can write 


r H (k) = r 0 {k) + 6(k) cosa, (50) 

where r Q and 6 are ordinary functions of the magnitude of the neutrino momen¬ 
tum and a is the angle between k and B. For a certain value of k, Ecp m 
determines a spherical shell limited by the spheres of radii r 0 {k) ± S(k), where 
the resonance condition is verified for neutrinos (or antineutrinos) moving in 
different directions with respect to the magnetic field. The quantity r 0 corre¬ 
sponds to the radius of the resonance sphere when B = 0 and in the case of 
oscillations between v e and v p (or v T ) is given by 

GfV2Y £ -^— = ± cos 2 6 . (51) 

m n To 2k 

The left-hand side corresponds to the difference b e — b pT evaluated at r = r Q and 
is positive definite. Therefore, assuming Am 2 cos 29 > 0, the above condition 
can be satisfied for u e <—> v p T but not for v e <—> To be effective as a 

kick mechanism the resonance transformations of active neutrinos have to take 
place in the region between the electron and muon (tau) neutrinospheres. In 
this case, the v e are trapped by the medium, but the v p j produced this way are 


1/3 
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above their neutrinosphere and, when moving outside, can freely escape from 
the star. In the presence of a magnetic field, the emission points for having 
the same k but different directions are not at the same radius, which originates 
an asymmetry in the momentum they carried away. 

Replacing k in Eq. <eh> by the thermal average (k) — 3.15 T and taking into 
account that p ~ 10 11 gem 3 and T ~ 4MeV at _R„ e , the requirement of having 
the resonance within the electron neutrinosphere implies Am 2 cos 28 > 2 x 
10 4 eV 2 , which is excluded by the experimental results on neutrino oscillations 
and the cosmological limits. For this reason, in what follows we concentrate 
only on the active-sterile oscillations. In this case, the resonance condition 
when B = 0 takes the form 


k(r 0e ) = T-^r-cos26>, (52) 

for <—> v s (ye <—* v s ). In the static model we use for the protoneutron star, 
a good approximation is to take Y e < 1/3. Therefore, if Am 2 cos 28 > 0, the 
resonant condition is verified only by antineutrinos and from now on we restrict 
ourselves to this situation. 

To determine the quantity 5 in Eq. ifolill we substitute Ve as given by Eq. 
in Eq. (EH) and expand be up to first order in 5e- Proceeding in this way we get 


with 


Se(k) = T>t(k) eB , 


W) = 7 ” 

K be 


(53) 

(54) 


where we have defined h g 1 = ^ hig for any function g(r). Using Eqs. ( 141111 - 1491 
we find explicitly 

V, = - -- 


2 


l Y. 


V„ tT = 


(1 Y e ) h p Y e h 


Ye 


/3 + Y e 

m, n Pe 

\ 

\5 Tm n 

4ir 2 p 

J 

3 T Y e 

m n p, e \ 


5 Tm n 

47 T 2 P ) 

r ° 


(55) 

(56) 


Finally from Eq. EEJ>, discarding terms proportional to the magnetic field, 
the following restriction to the vacuum mixing angle results 


tan "28 > 


b e 


be 


(57) 


which guarantees that the evolution will be adiabatic and an almost complete 
flavor transformation for small values of 8. 
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5 Neutrino momentum asymmetry 

If the temperature and density profiles are isotropic, then the momentum flux in 
the protoneutron star atmosphere is radial (see Eq. (|25l)). As a consequence, in 
most of the works on the kick mechanism driven by matter neutrino oscillations 
the resonance condition was evaluated assuming that neutrinos move outside 
with a momentum ( k) = (77r 4 /180 £(3)) T ~ 3.15 T pointing in the radial di¬ 
rection. In presence of a magnetic field this approach leads to the concept of a 
deformed resonance surface, which acts as an effective neutrinosphere. However, 
at each point within the neutrinosphere there are neutrinos with momentum k 
pointing in every directions and the resonance condition as given by Eq. 133 
actually defines a spherical surface for each value of cos a. This surface acts as 
a source of sterile neutrinos produced through the resonant Dp transformations, 
which move in the k direction. Therefore, as was mentioned above, for a given 
k the resonant transformations take place within a spherical shell, and not on a 
deformed spherical surface. A more careful calculation of the pulsar kick along 
this line has been carried out in Ref. and here we follow the same approach. 

Let us consider the situation where the resonance shell is totally contained 
within the Dp neutrinosphere. The sterile neutrinos produced there that move 
out freely escape from the protostar, but those directed toward the inside cross 
the resonance surface again and reconvert into Dp which, being within their own 
neutrinosphere, are thermalized. Consequently, only those D s going outward 
can leave the star and the resonance surfaces for different momentum directions 
behave as effective emission semispheres of sterile antineutrinos. For k pointing 
in opposite directions the radii of the respective semispheres are = r 0 (k) ± 
6(k) |cosa|, and therefore they have different areas that generate a difference in 
the momentum carried away by the sterile neutrinos leaving the star in opposite 
directions. There also exists a temperature variation within the resonance shell 
such that the temperatures of the emission semispheres for D s going outside in 
opposite directions are different (T(r“) > T(r+)). The variations in the area 
and the temperature tend to compensate each other when the contributions of 
D s emitted in every direction are added up. In the calculation below we take 
both semispheres into account when explicitly computing the asymmetry in the 
total momentum /C emitted by the cooling protoneutron star. In this way, we 
improve the calculation of Ref. m where the temperature within the resonance 
shell was assumed to be uniform. 

The active neutrinos and antineutrinos are emitted isotropically from their 
respective neutrinospheres and do not contribute to the momentum asymmetry 
A/C. Therefore, the only nonvanishing contribution can come from the D s . If 
their emission lasts for an interval At of the order of a few seconds, then in 
the static model we are using for the protostar A AC = K B At, where K B is 
the component along the magnetic field direction of the neutrino momentum 
emitted per time unit. According to Eq. m we can write 



( 58 ) 
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where k • B = cos 9 cos 9' — sin 9 sin 9' sin <f>' and the distribution function for 
sterile neutrinos f Bg is evaluated at r R . We have chosen a reference frame (x, 
y 1 z) fixed to the protostar where the magnetic field coincides with the z-axis, 
B = B z. In addition, at each point within the resonance shell we use the 
local frame introduced in Section 2 to evaluate the components of the energy- 
momentum tensor. Here, the upper limit for the integration on 9' is 7r/2 instead 
of 7 r, because we have to include only the contributions from the sterile neutrinos 
going outside. 

To proceed further with the calculation of K B we need to know / . The 
proper treatment of the interplay of collisions and oscillations for trapped neu¬ 
trinos requires the use of the density matrix. For our purpose a less elaborate 
description will suffice and we put = V{yt —> P s )/ P , where V(Pt —> D s ) is 
the probability for D( to convert into D s . Let us assume that the adiabaticity 
condition given by Eq. <E3 is verified in the interior of the vt neutrinosphere, 
and hence for small mixing V(yi —> v s ) = 1. Accordingly in Eq. (I58H we take 

/*>«) - /Jb). with 


f B .{r R ) = ftHr ol ) - 


cos 9' 


df eq 

"a 


ye. k 2 p dr 


df eq > 

-k-B -p-6 t (k) 
ar 


(59) 


where a term proportional to the derivative of p 1 df f : q /dr has been neglected, 
in agreement with the diffusive approximation we are employing. Substituting 
the expansion into Eq. and keeping terms that are at most linear in 
5e we find 


2 

K b =— / dk k 3 r oe (k) 6e(k) 

( f: q - 


df; q 


2x- k 2 p dr 


1 df eq > 
-r oe (k)-^6 e (k) 


(60) 


The energy released during the core collapse of the progenitor star is approxi¬ 
mately equipartitioned among all the neutrino and antineutrino flavors. In the 
context of our analysis this means that K r At ~ /C/6, where I\ r represents the 
momentum rate of the sterile neutrinos in the radial direction and is given by 
an expression similar to Eq. with k ■ B replaced by k • r . Following the 
same procedure outlined above we find 




df eq ' 


pk 2 dr 


(61) 


In terms of K B and K r the fractional asymmetry in the total momentum 
neutrinos carry away becomes 


A JC = \Kb\ 
1C 6 K r ' 


(62) 
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To evaluate the remaining integrals in Eqs. m and m the explicit dependence 
on k of the functions r oc and T>z have to be known. Simple analytical results 
can be derived by replacing these functions by the constant quantity f 0e and 
T>i = T>i(f oe ), where f oe is determined from the resonance condition when k is 
replaced by its thermal average 3.15 T(r): 

6.3 T(f ol ) bi{f ot ) = Am 2 cos 29. (63) 


Proceeding in this manner we get 


K r = 

K b “ 



7tt 2 T 4 

4 dT 2 \ 

48 ^ 

5 

3k dr J 

nr 0t 5t 

( 7tt 2 T 4 

7iT 2 f 0t dT 4 

36 

5 

10 dr 


1 dT 2 \ 
Xp e P dr ) 


(64) 

(65) 


In these expressions the terms that depend on dT 2 /dr come from the diffusive 
part of the distribution function. In the regime considered here they are smaller 
than the other contributions and in what follows we will neglect them. 

Suppose that the pulsar momentum K pu i is entirely due to the active-sterile 
oscillation driven kick, then A K, = K pu i. For a pulsar with a mass M « Mq 
and a translational velocity v ~ 500 kms -1 we have K pu iK, 10 41 gems -1 . Since 
neutrinos are relativistic, the total amount of momentum they carry is equal to 
the gravitational energy liberated by means of them, that is K, ss 10 43 gems -1 . 
Therefore, as mentioned in the Introduction, A/C//C must be of the order of 0.01 
to obtain the required kick. In this case Eq. El yields 


eB = 


0.045 f Q 


K 1 


- r o; hrp 




( 66 ) 


with h y 1 denoting the logarithmic derivative of T[r) evaluated at r 0e ■ This 
formula allows us to determine the magnitude of the magnetic field once the 
temperature, the baryon density, and the electron fraction profiles are known. 
Below we do it for the v p —v s conversion under the assumption that the resonance 
layer is entirely within the P M -sphere. However, we will first indicate how Y e (r) 
is determined by means of T(r) and p (r) within our approximate description of 
a protoneutron star. 

For p, ~ 0, the (3 equilibrium gives the relation p n — p v ~ p e among 
the chemical potentials of the background fermions. Hence, for non-relativistic 
nucleons N n /N p ~ exp[— (m n — m p — p e )/T] and 


Y, ~ 


1 


1 + e^“l T ’ 


(67) 


where we used the fact that in the condition prevailing in a protoneutron star p, e 
is much larger than the difference between the neutron and the proton masses. 
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Since electrons are degenerated p, e = (37T 2 pY e / rrinf ^ and Eq. is an im¬ 
plicit equation for Y e to be solved numerically once p and T are known as a 
function of r. Taking the derivative of Eq. EH with respect to r, h Y can be 
expressed in terms of h~ x and h Y l as follows 

h Y l = - ( h p 1 - 3 hj, 1 ) ^ _y _)_ 3 t ■ (68) 

Me 

To estimate B we now assume that the density and the temperature are 
related according to p = p c (T/T c ) 3 where p c and T c are the respective 
values of these quantities at the core radius R c . From this expression we see 
that h^ 1 = 3 h Y which implies that h Y ^ = 0 (see Eq. (El). Therefore, the 
electron fraction Y e and as a consequence the ratio /r e /T are constants. As core 
parameters we take R c = 10 km, p c = 10 14 gcm -3 , and T c = 40MeV, while for 
the density profile we adopt the potential law p(r) = p c (f? c /r) 4 which yields 
h^ 1 = —4 /r and T(r) = T c (R c /r ) 4 / 3 . The i^-sphere is taken stationary at a 
density of about 10 12 gcm -3 that corresponds to a radius ~ 3 R c and a 
temperature T R _ ~ 9MeV. Using these results as well as the formula for V R 
given in Eq. El , from Eq. El we obtain the following simple result: 

B = 5 x 10 17 (l + 0.3^ |-G, (69) 

where T = T{f 0 p ), with R c < f 0 < We see that the required intensity 

of the magnetic field decreases monotonically from 3.8 x 10 1 ' G at the core 
radius to 4.8 x 10 16 G at the ^-sphere. According to Eqs. El and im . within 
the specified interval for f 0 the range of the allowed values of the oscillation 
parameters are 7 x 10 s eV 2 > Am 2 cos 29 > 10 5 eV 2 , and for small mixing we 
have 25KeV > m s > 0.3 KeV. The requirement that the resonant conversion 
be adiabatic imposes the condition tan 2 28 > 3 x 10 -11 ( r 0l ,/Rc ) on the mixing 
angle. A last comment is in order, the ratio of the diffusive term to the isotropic 
one in Eqs. El and (l?)51i can be expressed as 10 4 (r 0(l /f? c ) 17 / 3 , which illustrates 
the validity of the approximation done in writing Eq. m . In the next section 
we repeat these calculations in the realm of the spherical Eddington model, 
which provides a simple but quite physically reasonable description of a static 
neutrino atmosphere. 


6 Spherical Eddington Model 

The actual configuration of the neutrinosphere and the resonance region de¬ 
pend on details of the protoneutron star, such as the density, temperature, 
pressure, and leptonic fraction profiles. The asymmetry in the momentum flux 
also depends on these features and on the configuration of the magnetic field. 
The present knowledge about protoneutron stars is not enough to single out a 
definite model for their structure, but there are some well established general 


17 



characteristics that help towards proposing models that are both workable and 
plausible. An important feature is that the size and the shape of the different 
profiles do not suffer sudden changes during the time interval when the main 
neutrino emission takes place, i.e. 0.5 s < t < 10 s. Thus, as a first approxi¬ 
mation we can consider a stationary regime that greatly simplifies calculations. 
Another characteristic to be taken into account is that the bulk of neutrinos 
is produced in the inner core. Consequently, we assume that there is no neu¬ 
trino production outside the core and hence the flux of diffusing neutrinos is 
conserved in this region. This flux is responsible for the transport of the energy 
liberated by the protostar. 

As mentioned before, the system is supposed to be constituted by nucleons, 
with a density p, electrons, neutrinos and antineutrinos. Neutrinos are in ther¬ 
mal equilibrium with nucleons, satisfying a transport regime consistent with the 
diffusion approximation, and are assumed to have a vanishing chemical poten¬ 
tial, /Tj, = pp = 0. The total energy density, the radial energy flux, and the 
pressure for neutrinos and antineutrinos in the case of spherical symmetry, can 
be obtained from Eqs. (El) . El) and mu 


u = J2 u „ e 

i 

i 


7 t 2 T 4 
-7T 1 , 

40 


1 dT 2 
36 xp dr 

U 

3"’ 


(70) 

(71) 

(72) 


where x 1 = x 1 + 2x 1 . The main contributions to the neutrino opacities 

ve i-'/a, T 

were discussed in Section III. Using the results given there for x ue and x v t 
(E qs. 1211) and 1221) ) we obtain x = 3 x 10~ 26 MeV -5 . 

If in addition we assume that baryons constitute a nonrelativistic ideal gas, 
we have the Eddington model. It is simple and physically well justified, al¬ 
lowing a detailed discussion of the relevant characteristics of the protostar and 
the geometry of the resonance region. Originally proposed to describe a stellar 
photosphere, this model was adapted by Schinder and Shapiro[IB] to a neutrino 
atmosphere with a plane geometry. The more realistic case of a spherical at¬ 
mosphere was considered in Ref. D1 to analyze the geometrical effect on the 
asymmetric neutrino emission induced by resonant oscillations. Photons and 
electrons are of course present and, in fact, electrons make the leading contribu¬ 
tion to the effective potential in the case of oscillations among active neutrinos. 
However, we can ignore both photons and electrons for the hydrodynamical 
description of the system. 

The state equation that defines the system is given by the sum of the con¬ 
tributions from nucleons and neutrinos to the pressure gas, which is 


m n 3 


(73) 
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In the Newtonian limit for the metric ( g oa = — 1 — 2 (f>, g a j = 0, and (jij = Sij , 
with (f> = — GM (r)/r) the conservation of the total (neutrinos + matter) energy- 
momentum tensor yields the energy flux conservation 


d(r 2 F) 

dr 


= 0 


and the hydrostatic equilibrium equation 


(74) 


dP 

dr 


(.P + p+U) 


GM 


(75) 


Here, M(r) = 4-7T JJ dr' r' 2 p{r') is the mass enclosed up to a distance r from 
the center. 

Through the atmosphere the baryon density is p ~ (10 11 — 10 14 )gcm -3 
and T ~ (4 — 40) MeV, and thus we have U — (10 -3 — 10 -2 ) p and P ~ 
(4 x 10 -3 — 4 x 10 -2 ) p. This means that on the right-hand side of Eqs. (17311 and 
d we can ignore the contributions coming from the energy density of neutrinos 
and the baryonic density. In addition, Eq. d implies that F = L c /47rr 2 , where 
L c is the luminosity of the protostar. Taking all the above into account we 
arrive at the following set of equations for an isotropic atmosphere of neutrinos 
in thermal equilibrium with an ideal gas of nucleons: 


dT 2 9x L c 

, = P 2 ’ 

ar nr* 

(76) 

P = — T, 

(77) 

m n 

dP GM 

(78) 

dr ^ r 2 ’ 

dM 2 

—— = Airr p . 
dr 

(79) 


Solving this system of equations yields the profiles of four functions through¬ 
out the atmosphere: pressure P(r), temperature T(r), baryonic density p(r), 
and the enclosed mass M(r). From them any other functions of interest can 
be calculated. To analyze the system it is convenient to introduce adirnen- 
sional variables normalized to the corresponding values at the core: x = r/R c , 
m(x) = M/M c , t{x) = T/T c , g(x) = p/p c , and p(x) = P/P c - Proceeding in this 
way, Eqs. d-d can be rewritten 


dt e 

dx °tx 2 ’ 

(80) 

p = tg, 

(81) 

dp mg 

~1 C C 2 ” , 

ax x* 

(82) 

dm 2 

— = d c x Q, 
dx 

(83) 
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where 


b r = 


9 xp c L c 
2nR c T? ’ 


Gm n M c 

RcT c 


d r = 


AiTp c Rl 

M r _ 


(84) 


are adimensional constants. Thus the set of functions that solve the previous 
system of equations depends only on these three combinations of the core pa¬ 
rameters. 

By eliminating g in Eqs. and m we arrive to 


dp c c dt 
dx b c dx 

The last equation is easily integrated as follows: 


P = 


t 2 — a 
1 — a 


where 


i{x) = 1 — 


2 b c . 


c c m 


( x ) 


and m(x ) is an effective mass defined by 

r, , , ,, dt _,./■* dt 

/ dx m(x ) t —— = m(x) / dx t —— 
J i dx J i dx 


(85) 

( 86 ) 

(87) 

( 88 ) 


From Eqs. m and we can also express the density in terms of the tem¬ 
perature as 

t 2 — a 

Q = f(l — a) ’ (89) 

which, together with Eq. m , yields a first order differential equation for t 

dt b c t 2 — a 


dx t 2 x 2 1 — 


= 0 . 


(90) 


At x = 1 (R = R c ) we have dt/dx | x=1 = — b c , which is independent of a (a;). 
For an idealized atmosphere, where the model would apply to the whole space, 
the physical solutions would correspond to an infinite protostar where the tem¬ 
perature has an asymptotic behavior for x 1, such that the adimensional 
temperature tends to t s ~ \fa. Thus, the function a(x) varies in the range 
1 — 2 b c /c c < a < t 2 for 1 < x < oo. 

Eq. 19011 has no analytical solution when a is a function of x. An approximate 
solution can be found by the following procedure. Let us consider the differential 
equation m with a constant. Then, an analytical (implicit) solution is given 
by 


t ~ 1 + ~2 


In 


— In 


]__ya\ 
1 + / 



(91) 


with 1 > t > \/a. If we replace the constant a by a well behaved function of x, 
the above expression still satisfies Eq. pi at x = 1. For an infinite atmosphere, 
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a good approximation to the exact solution is given by Eq. ED, with a now a 
function of t(x): 

a(t) = 1 - — - A(1 -t) , (92) 

c c 

where A = ^1 — / (1 — t s ). From Eq. (19211 . we see that a(t s ) = t 2 s and 

a(l) = 1 — 2& c /c c , which means that this ansatz fits the extreme values of a(x). 

Once p and T are known, we can find the electron fraction Y e as a function 
of r by solving the implicit expression derived in Eq. ED- To examine the 
prediction of the model in connection with the kick mechanism, we consider 
a protoneutron star with reasonable values for its core parameters, the same 
as in Ref. '3- Thus we take M c = = 1.989 x 10 30 kg, R c = 10 km, 

L c = 2 x 10 52 ergs _1 , p c = 10 14 gcm -3 , and T c = 40MeV. In this case, 

b c = 1.96, c c = 3.52, d c =0.625. (93) 

A numerical estimation for the asymptotic value of the solution of Eqs. ED-EP 
gives t s = 0.12. The electron fraction takes the value Y e = 0.08 at the core. 

According to Eq. rcni in terms of the adimensional functions the mean free 
path for the electron neutrino becomes 


g 

=^^cm. (94) 

Inserting the above expression in the integrand of Eq. m and evaluating the 
resulting integral numerically we obtain x Ve ~ 2.6 for the radius of the v e - 
sphere. This yields p(R Ve ) ~ 1.3 x 10 n gcm -3 , T(R Ue ) ~ T s = 4.8MeV, and 
Y e (R Ve ) ~ 0.1, which are in reasonable agreement with typical values from 
numerical simulations l?21 ffl . In the same way, for the muon (tau) neutrino we 
get x v T ~ 2.2. The profiles that correspond to the solution of this model are 
given in Fig. 1. 

From Eq. ED we can estimate again the intensity of the magnetic field re¬ 
quired for the kick in the case of 4 <—> D s . Besides the temperature, density, 
and electron fraction profiles we also need their logarithmic derivatives. For T 
and p they are immediately calculated from Eqs. ED-ED, and the correspond¬ 
ing results are 


Ziy 1 = 


- 6 , 


h~ l + h? 1 = 


Q_]_ 

Ax r ’ 

m 1 

tx r 


(95) 

(96) 


The formula for h Y l is obtained substituting the above expressions into Eq. ED- 
Now we have all the quantities needed to determine the magnetic field from 
Eq. IH?)1) in the spherical Eddington model. The result is plotted in Fig. 2. 
From the curve presented in this figure we see that B decreases from a value of 
8 x 10 16 G in the core up to a minimum of 3.8 x 10 16 G at 7" 0 ~ 16 km. From 
this point the magnetic field increases steadily as we approach the surface of 
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the neutrinosphere. In contrast with the model used at the end of Section 5, 
here at r ~ 21 km the contributions from the geometrical and the temperature 
variations compensate each other. At this radius the factor (1 + 2f 0t ) in the 
denominator of Eq. m vanishes and the magnetic field takes arbitrarily large 
values. In other words, no kick can be generated by this mechanism in the 
regions near the surface of the i/^-sphere. 

In the region between the core and the surface of the neutrinosphere R c < 
r Q < R v , the range of values allowed for the oscillation parameters is 

10 9 eV 2 > Am 2 cos 2 6 > 10 6 eV 2 , (97) 

For a small mixing angle this inequality translates into the condition 30 KeV > 
m s > 1 KeV for the sterile neutrino mass. Small mixing is allowed by the 
adiabaticity condition, which warrants an efficient flavor neutrino conversion 
for tan 2 29 > 7 x 10~ 12 at r 0(l = R c and tan 2 29 > 3 x 10 -9 at r 0(l = R V/J . The 
results for m s and the mixing angle agree with those derived in the previous 
section and are compatible with a sterile neutrino that is a viable dark matter 
candidate@2j. 

7 Conclusion 

In this work we have examined in detail a possible explanation of the large drift 
velocities observed in pulsars in terms of an asymmetric neutrino emission. The 
asymmetry is a consequence of the resonant neutrino conversions affected by 
the strong magnetic field characteristic of protoneutron stars. We have shown 
that for active-sterile neutrino oscillations this is a feasible kick mechanism. 
Two conditions must be simultaneously satisfied to generate a natal kick by 
this mechanism. First, the conversion has to take place at the interior of the 
neutrinosphere of the active neutrino, and second the magnetic field has to 
be intense enough in order to induce a momentum asymmetry of the required 
magnitude. 

We have made detailed calculations by means of the neutrino distribution 
function in the diffusion approximation, combined with the idea of a spherical 
resonance shell for neutrinos with momentum k that move in different directions 
relative to B. This scheme provides a better description of the problem than the 
one formulated in terms of a single deformed sphere acting as an effective emis¬ 
sion surface of sterile neutrinos. In our approach the sterile neutrinos produced 
in the resonance shell move in every direction relative to B. However, those 
going toward the interior of the protoneutron star cross again the resonance 
region and are reconverted into active neutrinos. These are within their own 
neutrinosphere and become thermalized. As a consequence, only the outgoing 
sterile neutrinos contribute to the total momentum emitted by the protostar. 

Two opposite effects have to be taken into account when computing the 
fractional momentum asymmetry. One is purely geometric, and comes from the 
difference in the areas of the semispherical emission surface for neutrinos moving 
in opposite directions with the same momentum magnitude. The other is due to 
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the radial temperature gradient, by which sterile neutrinos produced at different 
depths have unequal energy. In general both contributions do not compensate 
and thus we obtain a non null fractional asymmetry in the total momentum. 
Explicit results have been obtained for two simple models of a protostar, with 
the matter background assumed to be composed by nonrelativistic nucleons 
and degenerated electrons. One of them, which we worked with in detail, is the 
Eddington model for a spherical neutrino atmosphere. 

Magnetic fields of the order of 10 16 —10 17 G are needed to reproduce the 
observed pulsar velocities. At first sight these fields could seem rather large, 
compared with the estimated values at the surface of a protoneutron star ( B s ~ 
10 13 G), but in fact intensities as high as B c ~ 10 18 G are possible at the core. 
A given parametrization for the magnetic field at the interior of a protoneutron 
star is|3B] 


B — B s + B c 


l _ e —/3(P/P,) 7 


(98) 


with (3 ~ 10 -5 and 7 ~ 3. By adopting this profile as an upper boundary for 
the magnetic field in the protostar, we can see that the required field remains 
well below this boundary in most of the region within the neutrinosphere. In 
addition, the oscillation parameters are compatible with the allowed region for 
sterile neutrinos to be warm dark matter, leaving the mechanism as an attractive 
possibility to explain the proper motion of pulsars. 
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Figure captions 

FIG. 1: The different profiles normalized to the corresponding value at the core 
for the Eddington model discussed in Section VI. 

FIG. 2: The magnetic field B required to produce the kick by the —> v s 
resonant conversion for the Eddington model discussed in Section VI, normalized 
to the field required at the core, B c . 
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